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Q\ • Abstract 

For the unitary ensembles of N x N Hermitian matrices associated with a weight 
function w there is a kernel, expressible in terms of the polynomials orthogonal with 
respect to the weight function, which plays an important role. For the orthogonal and 
symplectic ensembles of Hermitian matrices there are 2x2 matrix kernels, usually 
constructed using skew-orthogonal polynomials, which play an analogous role. These 
matrix kernels are determined by their upper left-hand entries. We derive formulas 
expressing these entries in terms of the scalar kernel for the corresponding unitary 
ensembles. We also show that whenever w'/w is a rational function the entries are 
equal to the scalar kernel plus some extra terms whose number equals the order of 
w'/w. General formulas are obtained for these extra terms. We do not use skew- 
orthogonal polynomials in the derivations. 

0\ 



1. Introduction 

o ' 

In the most common ensembles of N x iV Hermitian matrices the probability density 
Pn(xi, ■ ■ ■ , %n) that the eigenvalues lie in infinitesimal neighborhoods of X\, ■ ■ ■ , Xn is given 
by 

P N (X!, - ■ ■ ,X N ) = C N Y[ \Xj -Xkf YlwiXj), 

j<k j 

where (3 = 1, 2 or 4 (corresponding to the orthogonal, unitary and symplectic ensembles, 
respectively), w(x) is a weight function and cjv is a normalization constant. 

For the unitary matrix ensembles an important role is played by the kernel 

JV-l 

KN(x,y)=^ f (pk(x)(p k (y), (1.1) 

k=0 

where {cpk(x)} is the sequence obtined by ort honor malizing the sequence {x h w(x) 1 ^ 2 }. The 
probability density is expressed in terms of it by 

P N (xi, ■ ■ ■ , x N ) = — det (K N (xj,x k ))j,k=i,~,N- 
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More generally the n-point correlation function R n (xi, • • • , x n ), the probability density that 
n of the eigenvalues, irrespective of order, lie in infinitesimal neighborhoods of x\, • • • , x n , is 
given by the formula 

Rn{xi, ■ ■ ■ ,x n ) = det {K N (xj, x k ))j,k=i,-, n - (I- 2 ) 

And the probability E(0; J) that the set J contains no eigenvalues is equal to the Fredholm 
determinant of the kernel K^(x,y) Xj(y), where X denotes characteristic function. 

For the orthogonal and symplectic ensembles there are 2x2 matrix kernels which play 
analogous roles. In this case the determinant in (|1.2|) is to be interpreted as a quaternion 
determinant (it is a linear combination of traces of products of matrix entries of the block 
matrix on the right side), and the square of E(0; J) equals the Fredholm determinant of the 
matrix kernel. (The last fact can be deduced from the computation in ||, sec. A. 7. A direct 
derivation is given in ||.) In the case of the orthogonal ensembles we shall always assume 
that N is even. The kernels for the orthogonal and symplectic ensembles are of the form 

/ Sm(x,y) S N1 D(x,y) \ 

K m (x,y)=\ (1.3) 
\ISm(x,y)-e(x-y) S N1 (y,x) ) 

and 

, / S m (x,y) S N4 D(x,y) \ 

K m (x,y) = -\ (1.4) 

\ IS m (x,y) S N ^(y,x) J 

respectively. Here e(x) = \ sgn(:r) and the explanation for the notation is this: the Snp(x, y) 
are certain sums of products and if Sn/3 is the operator with kernel Snp(x, y) then Sn/3D(x, y) 
is the kernel of S^pD (D = differentiation) and IS^pix^y) is the kernel of I Sn/3 (I = 
integration, more or less). We shall write these out below. One can see from this description 
that once the kernels Sjsr^x^y) are known then so are the others. 

Matrix kernels were first introduced by Dyson 0] for his circular ensembles and he es- 
tablished the analogue of formula ( |1.2| ) for the correlation functions. Later, Mehta and 
Mehta and Mahoux found matrix kernels for the ensembles of Hermitian matrices, and ex- 
pressed them in terms of systems of skew-orthogonsal polynomials. These are like orthogonal 
polynomials but the inner product (different in the (3 = 1 and (3 = 4 cases) is antisymmetric 



instead of symmetric. In terms of them one obtains for S^pix, y) sums like the one in (|1.1| ) 
but which are a little more complicated. A problem here is that the skew-orthogonal poly- 
nomials are not always that easy to compute and, even if they are, the sums involving them 
may not be easy to handle. For example, one is often interested in scaling limits as iV —>■ oo 
and in order to do this it helps to have a good representation for the sum. 

In this paper we shall not use skew-orthogonal polynomials at all. Instead, we shall use 
the expressions for the various matrix kernels in the general form given in M, and derive 
general formulas for the Sn/3{x, y) in terms of the scalar kernel K^{x, y) given by ( |1 . 1|) , with 
N replaced by 2N when (3 = 4 and w replaced by w 2 when (3 = 1. More exactly, we shall 
express the operators whose kernels are the S^pix, y) in terms of the operator whose kernel 
is Kjy(x,y). These are given in Theorem 1 below. 
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The formulas can be brought to a very concrete form whenever the support T> of w 
is a finite union of finite or infinite intervals and w'/w is equal to a rational function on 
T>. (Such weight functions are called semi-classical since they include the weight functions 
for all the classical orthogonal polynomials.) We find then that the S^p{x,y) are equal to 
the appropriate scalar kernel Kj^{x : y) plus some extra terms whose number is independent 
of N. This number equals the order of w'/w, the sum of the orders of its poles in the 
extended complex plane. We must also count as a simple pole any end-point of T> where 
w'/w is analytic. Thus, for the Gaussian ensembles (w(x) = e~ x ) and Laguerre ensembles 
(w(x) = x a e~ x ) there will one extra term because of the simple poles at oo and respectively, 
and for the Jacobi ensemble (w(x) = (1 — x) a (l +x)") there will be two extra terms because 
of the simple poles at ±1. For the Legendre ensemble on (—1, 1) there will also be two extra 
terms although w'/w = in this case. 

We shall produce explicit formulas for the extra terms, which are given in Theorem 2. 
These will be used to work out the cases of the Gaussian ensembles (well-known 0]) and the 
Laguerre ensembles (known apparently only in the case a = ||). 

To apply our formulas to the Laguerre ensemble we require at first that a > so that 
Theorem 1 is applicable. The formulas for general a > — 1 are then obtained by analytic 
continuation. Similar analytic continuation arguments apply quite generally. (See the remark 
at the end of section 3.) For example, for the Legendre ensemble we would start with the 
formulas for the Jacobi ensemble for a, (3 > and then take the analytic continuation (or 
limit) to obtain the fomulas for a — [3 — 0. This is the reason the end-points ±1 count as 
poles. 

The recent announcement || has some elements in common with ours. A generalization of 
the Laguerre ensemble was considered there where e~ x was replaced by the exponential of an 
arbitrary polynomial and the occurrence of only finitely many extra terms was established, 
without their being evaluated, using skew-orthogonal polynomials. This fact was used to 
deduce universality for this class of ensembles. 



2. The general identities 

We start with the expressions for the various matrix kernels in the form given in || . (The 
notation here is slightly different.) Taking the symplectic enembles first, we let {pj(x)} be 
any sequence of polynomials of exact degree j and define (fj(x) = Pj(x) w(x) 1 ^ 2 . Let M be 
the 2N x 2N matrix with j, k entry (j, k = 0, • • • , 2N - 1) 

m jk = ^J {Pj{x)p' k {x) -Pj(x)p k (x)) w{x) dx = ^J (<pj(x) y' k {x) - tp'jix) (p k (x)) dx. (2.1) 

This matrix is invertible and we write M~ l = {fijk)- Then 

2N-1 

S Ni (x 1 y)= <Pj( x )f J 'jk<Pk(v) (2-2) 

j,k=0 

and 

2AT-1 27V-1 

IS N4 (x,y)= Vjk<Pk(y), S NA D(x,y) = - ^ (pfa) fi jk <fi' k (y). 

j,k=0 j,k=0 
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Any family of polynomials leads to the same matrix kernel. Of course at this point the 
formulas look quite bad because of the 

For the orthogonal ensembles we take we take the pj as before but this time define 
(fj(x) = Pj(x) w(x) and let M be the N x N matrix with j, k entry (i, j = 0, • • ■ , N — 1) 



m jk = J j e(x -y)p j (x)p k (y)w(x)w(y)dydx = J ipj(x) eip k (x)dx. (2.3) 

Here e denotes the operator with kernel e(x — y). Again M is invertible, we write M _1 = 
(yUjfc), and the formulas for the kernels are 

N-l 

S m (x,y) = - <Pj( x ) Hk£<Pk(y), 

j,k=0 

N-l N-l 

IS N1 (x,y) = - £( Pj( x ) Vjk£<Pk(y), S N1 D(x,y) = J2 <Pi( x ) A*J* <pM- 

j,k=0 j,k=0 

We change notation so that we can treat the two cases at the same time — we shall see 
that they are interrelated. We continue to use the notations N and w, but when {3 = 4 the 
iV here will be the 2N of ( |2.2j ) and when (3 = 1 the w here will be square of the weight 
function in (|2.3|) . Thus in both cases iV is even, we take pj to be polynomials of exact degree 
j and set tpj = pjw 1 ^ 2 . The matrices {w!f^) and are the M and M -1 corresponding 

to the (3 — 4 and 1 ensembles. We set 

N-l N-l 

$n \x,y) = £ ij$ (Pk(y), s ( j}\x,y) = - J2 vA x )Vjk £ Vk{y)- 

j,k=0 j,k=0 

Finally, K N (x,y) will denote the (3 = 2 scalar kernel (|1.1|) . 

We denote by TC be the linear space spanned by the functions (fo,- ■■jtpN-i, in other 
words the set of all functions of the form u> 1//2 times a polynomial of degree less than N. We 
denote by K be the projection operator onto TC. Its kernel is K^(x,y). Finally, we denote 
by the operator with kernel S$(x,y) and by S^' the operator with kernel Sj^\y,x). 

The following lemma will identify these operators. We think of our weight functions as 
defined on all of R, and our basic assumptions are 

HCL 1 {B), DHcL l (R). (2.4) 

The former is needed even to define the ensembles. The latter is restrictive and implies in 
particular that all the ip k are absolutely continuous. We use the notations D« and En for 
the restrictions of the operators D and e, respectively, to TC. 

Lemma. The operators KD^ and Ke-h are invertible and 

SW\ n = D{KD H )-\ S^\ H ,=0, 

S^\n=e(KE H )-\ S^V = 0. 
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Proof. Integrating by parts the second integral in (|2.1|) shows that mk = / <Pj(x) <p' k ( x ) dx. 
(This is where the absolute continuity of the tpk come in.) Thus for i — 0, • • • , N — 1, 

3,k j 

Since the ^ span H we see that K <p' = <p' for all ip E Tt. This shows that KD-^ is a 
one-one, and hence invertible, operator on H, and also that S^\n = D(KDu)~ l ■ Of course 
S^l-^x = since each e 7Y. This proves the first part of the lemma. For the second, 
observe that by the antisymmetry of (m^) the formula for S$(y,x) can be obtained from 

the formula for Sff(x,y) by replacing m$ by and <p'j(x) by Eipj(x). Thus the second 
part of the lemma can be proved just as the first, replacing D everywhere by e. 

To identify (KD-^) -1 and (Ken)' 1 more concretely we shall enlarge the domains of D 
and e. We have DH C -^ 1 (R) by assumption, and eH C L°°(R) since H C L 1 (R). It is easy 
to see that the operators 

D:H + eH^H + DH, e:H + DH^H + eH 

are mutual inverses. In the following, I-h+dh an d Ih+eH wm denote the identity operators 
on the spaces H + DH and H + eH, respectively. 

Theorem 1. We have 

S (4) = (Ih+dh -(I- K)DKe)- l K, (2.5) 

S {1) ' = (I n+ en -(I- K)eKD)~ 1 K. (2.6) 

Proof. Since D and e are mutual inverses we might guess that a good approximation to the 
inverse of KD-^ is Ke-n, where e-h denotes the restriction of e to H. With this in view, we 
compute 

KeKD h = KeD h -Ke(I- K)D H = I n -Ke{I- K)D n , 

where denotes the identity operator on H. The operator on the right side is invertible 
since both KD n and Ke h are, and we deduce that 

(KDn)- 1 = (I H - Ke(I - K)D h )- 1 Ke h . 

Hence by the lemma, 

S^\ n = Dh(KDh)- 1 = KDn(KDn)- 1 + (/ - K)D n (KD H )- 1 

= I H + (I- K)D h (Ih - Ke(I - K)D u Y x Ke u . 
Recall that the domain of e is H + DH and set 

A = (I - K)D H :H^H + DH, B = Ke :H + DH -> H. 

Then 

Ih+dh + (I ~ K)D H {I H - Ke(I - K)D h )- 1 Ke 
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is equal in this notation to Ih+dh + Mjn ~ BA) 1 B. This in turn equals (Ih+dh — AB) 1 . 
Hence restricting to H gives 



S^\ n = {I n+DH - {I - K)DKe) 



H 



Since S^\^± = K\ n ± = this gives Q2. 
interchanging the roles of D and e. 



and (|2.6| ) is obtained by an analogous argument, 



Remark. The identities of the theorem may be restated in the rather more complicated 
form 

= K + (I- K)DKe(I H+DH -(I- K)DKs)- l K, 
SO)' = K + (I - K)eKD{I H+eH -(I- K)eKD)- l K. 

The summands on the right may be thought of as corrections and we see that they will be 
of finite rank (independent of N) whevever (J — K)DKe and (I — K)eKD are. This will be 
true whenever the commutator [D, K] is, which will be the case in what follows. 

3. The case of rational w'/w 

We assume now that w'/w is a rational function on the support of w and, at first, that 
Q2.4|) holds so that Theorem 1 is applicable. We explain at the end of this section how to 
remove the restriction in the cases of greatest interest. From now on it will be convenient 
to take the pj to be the polynomials orthonormal with respect to the weight function w so 
that the (pj are orthonormal with respect to Lebesgue measure. 

It follows from the Christoffel-Darboux formula that there is a representation 



K N (x,y) = a N 



<p N (x) <p N -i(y) - <Pn-i(x) (Pn(v) 
x-y 



/ 1 \ / <p N (y) 
= a N ((p N (x) <p N -i{x) ) j \/{x-y) (3.1) 

V -1 J \ cp^y) 

for a certain constant a^. This holds for an arbitrary weight function. Whenever w'/w is a 
rational function there is a differentiation formula 



A B 



-C -A 



Vn-i 



where A(x), B(x) and C(x) are rational functions whose poles are among those of w'/w, 
counting multiplicity. (See 0, sec. 6.) From this and ( |3.1| ) we find that the kernel of [D, K], 
which equals (d x + d y ) K^{x, y), is equal to 



( C(x) - C{y) A{x) - A(y) \ 
x — y x — y 

A(x)-A(y) B(x)-B(y) 




(3.2) 



x-y 



x-y 
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It follows from this that [D, K] is a finite rank operator, and that its kernel is expressible in 
terms of the functions 

x k ip N -i(x), x k ip N (x), (0 < k < Hoc) (3.3) 

where is the order of w'/w at infinity and, for each finite pole Xi of w'/w, the functions 

(x -Xi)" fc ~Viv-i(a;), (x -Xi)-*' 1 y N {x), (0<k<n Xi ). (3.4) 

where n x . is the order of w'/w at This is seen by expanding the functions appearing in 
the central matrix in (|3.2|) , which is simple algebra. 

In the space spanned by these 2n functions (n is the total order of w'/w) there is a 
subspace of dimension n contained in 7i and a subspace of dimension n contained in Ti.- 1 . To 
see the first, an inductive argument using the three-term recurrence formula shows that the 
subspace spanned by the funcions ( |3.3|) contains the n M functions (pN-k (0 < k < n^) which 
lie in Ti. The functions ( |3.4j ) span a space of dimension 2 J2 Ti Xi consisting of functions which 
equal w 1 ^ 2 times rational functions which may have poles at the Xi of order n Xi . A function 
in this space will belong to 7i if the principal parts at all these poles vanish. This gives X) n Xi 
conditions in a space of dimension 2 n xi , giving us a subspace of dimension X n Xi which is 
contained in TC. Thus the space spanned by the functions ( pT3| ) and ( |3.4j ) together contains 
a subspace of dimension n contained in H. To see that there is an n-dimensional subspace 
lying entirely in Ji , observe that 7i is spanned by the functions 

<PN-k {k < rioo), if k (k < n-Uoo), Yl(x - Xi) nx * x k (k < N - n). 

Our In functions are all orthogonal to the last of these, whereas orthogonality to the re- 
maining ones imposes n conditions, giving a subspace of dimension n which is contained in 

It follows from the preceding discussion that the space spanned by the functions (|3.3| ) 
and ( p.4|) contains n linearly independent functions ipi, ■ ■ ■ ,ip n lying in T~C and n linearly 
independent functions ipn+i, ■ ■ • , ip2n lying in Ti 1 . And we have a representation 

2n 

[D, K] = J2 Aj^i^^j (3.5) 

for some constants Aj which can be determined from ( |3.2| ) once we have fixed the (We 
use the notation a ® b for the operator with kernel a(x) b(y).) It follows that also 

2n 

[e, K]=J2 AjeA®^- (3-6) 

Here we used eD = Ds = I, the antisymmetry of e and the easy fact that (a® b) T = a®(T'b) 
for any operator T . These will be used again below without comment. 

The matrix A = (Aj) is symmetric since K is symmetric and D is antisymmetric. (We 
hope this A will not be confused with the function A appearing in ( |3.2| ).) Since K is the 
projection operator onto H the commutator [D, K] takes H to Ti 1 - and H 1 - to 0. Hence 

Aij = if i, j < n or i, j > n. (3.7) 
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After a little more notation we shall be able to state the formulas. We already have the 
matrix A. We define the matrix B by 

Bij = (eipi, ipj). 

Define J to be the matrix whoese i,j entry equals 1 if i — j < n and otherwise. Finally, 

set 

C = J + BA 

and write A for the matrix obtained from A by deleting its last n columns, Cq for the 
matrix obtained from C by deleting its last n rows and Coo for the matrix obtained from C 
by deleting its last n rows and its last n columns. Observe that by ( |3.7| ) the first n rows of 
A are zero. 

Theorem 2. We have 

S%\x, y) = K N (x, y)- E (A C m l (?„)<; ^(x) e^(y), (3.8) 

i>n, j 

S$(x,y) = K N (x,y)- E [AC (I - BAO)-% ^(x) e^iy). (3.9) 

i<n, j 

Proof. Using ( |3.5| ) we find 

(/ - K)DKe = [D, K]Ke = (E A v V»i ® Ke = ~ A v ^ ® #j 

since -K^i = when j < n and = when j > n. Thus 

I -{I- K)DKe = 1+ A n ^ ® e ^r 

j<n, i 

Now if we have a finite rank operator Yl, a i®^i then 

(J + E a, ® fei)- 1 = / - E T ^ a * ® 6 * ( 3 - 10 ) 

id 

where T is the matrix with entries 

Tij = dij + (bi, a,j). 

In our case i, j <n and 

ai = ^A ki ^k, bi = eipi, 

k 

so 

Tij = Sij + E^i, Ak i = % + X! 
This equals (J + BA)ij = Cij and so we have shown 

(/-(/- K)DKe)- 1 = / - E (Coo)7i (E^fc ® 

i,j<n k 
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whence 

= {I - {I - K)DKe)- l K = K - ]T £ A u (C o0 )7/^ ® Ke^. (3.11) 

i,j<n k 

To compute Keifij we apply ( |3.6|) to using the fact that x/jj G 7i, to obtain 
Ife^ = eif)j - A ik e^i (eif) k , 4>j) = #j - E A ik B kj e^i 

= #j + E( 5 ^W' = E CWi- (3-12) 

Here we used the symmetry of A and the antisymmetry of B. Substituting this into ( p.ll| ) 
gives 

S W = (/-(/- K)DKe)~ l K — K — £(A C" 1 C„)« V> fc ® e^,, 

it,i 

which is the same as Q3.8|). 

To derive (13.91) we use (13.61) and find that 



(/ - K)eKD = [e, K]KD = - £ ^ e^< ® Dtfe^-. (3.13) 

Using ( |3.5| ) again and the fact that De = I we see that 

DfTeVj = + E Vfc (Vk £, 0j) = + E ^ A ik^k- 

Again we use the fact that Kipj = ipj when j < n and i^V'i = when j > n. If we recall the 
definitions of J and C we see that we have shown DKeipj = J2kCjkipk- Substituting this 
into ( |3.13| ) gives 

(I - K)eKD = - ]T A ij C jk eipi ® ip k 

and so 

/ - (/ - K)eKD = 1 + J2( AC )*k eA ® Vfc- 
We use (p . 1 Q|) again, this time with i,j < 2n and 



Oi = ^(AC) fci eip k , h = A- 
k 

Now we have 

Tij = $ij + E(^> e^*) (^C)fcj = Sij - E ^ (-ACOfci = ( J " BAC)ij. 



Hence ( |3.10| ) gives 



(I-(I-K)eKD)- 1 = Y(I-BAC)^Y( A C)ki£^k®^ = ^[ACtf-BAd)- 1 ^^®^. 

i,j k j,k 
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To obtain S^' we must right-multiply by K, which has the effect of imposing the restriction 
j < n. After taking transposes and changing notation we obtain ( |3.9|) . 



Remark. Here is how to extend the results to the case where the second part of Q2.4j) may 
not be satisfied but the support T> of w consists of a finite union of intervals. Denote now by 
Xi the poles of w'/w together with all finite end-points of T> where w'/w is analytic. Then 
we can write 

w ( x ) = T[( x - x d ai w o( x ), 

i 

where wo satisfies ( |2.4| ) and each ai > — 1. Think of w, and therefore the kernels K^(x,y), 
S$(x,y) and S$(x,y), as functions of the a^. Theorem 2 would apply to w itself if all 
the «j > since then ( p.4|) would be satisfied. But the constituents of these kernels are 
real- analytic functions of the «j, so the formulas for aj > — 1 (and therefore for our given 
weight function w) can be obtained by analytic continuation of the formulas for > 0. 

4. The Gaussian and Laguerre ensembles 

These are (essentially the only) cases where n = 1 and are especially simple, as we shall 
now see. 

By the symmetry of A and (|3.7|) , A has the form 



.4 



/OA 



V A 



for some constant A. There arise two functions, tpi and fa, the first lying in 7i and the 
second lying in H L . 

Since the two 2x2 matrices A and B have diagonal entries, AB is a diagonal matrix 
and therefore so is C. Therefore ( |3.12 ), in which j = 1, says that Ksfa = Cufa. Since 



Ke^ is invertible C\\ ^ 0, and so efa G 7i. This implies that all entries of B vanish, so 
B = 0, C = J. It is immediate from these facts that 

/ 

A C 00 l C = AC(I-BAC)- 1 = 



V A 



Hence by Theorem 2, 



S ( N(x,y) = K N (x,y) - Xfa(x) £fa(y), 

S ( N(x,y) = K N (x,y) - Xfa(x) efa(y). 

It remains to find the constant A and the functions ip\ and fa in the two cases. 

The Gaussian ensembles. Here the pole is at x = oo. Clearly fa = (fiN-i and fa = <fiN 
in this case. Moreover we have for this ensemble 

a N = y/iV/2, A(x) = -x, B(x) = C{x) = V2N, 
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so we find from ( |3.2| ) that 



A 



\-Jn/2 



J 



which gives A 



N/2. Therefore 
s n\ x >v) = k n(x,v) + \jN/2 <p N (x) e(p N -!(y), 
S { N(x,y) = K N (x,y) + \Jn/2 (p N -i(x) e<p N (y). 



The Laguerre ensembles. Here w(x) = x a e x and 

Pj (x) = 



K.I n +l)^^' 



(4.1) 



where L,- is the generalized Laguerre polynomial. We assume at first that a > so that 



holds. In the notation of ( |3.1| ) and ( |3.2| ) we have 
a N = -y/N(N + a), A(x) = x~\N + |) - i B{x) = C(x) = -x^JnJn + a). 
The pole is at x = 0, and if we set £j(x) = x~ 1 ipj(x) then ( |3.2| ) becomes 



\J N(N + a) ( £ N (x) 6v-i(z) 



/ -yjN(N + a) 



N + 



a 



a 



V 



N + - —yN(N + a) 



\ 



J 



( Sn(v) 

V Civ-i(y) 



• (4-2) 



Our functions ipi and 4>2 are linear combinations of £tv and £at_i, with ipi lying in 7i 
and ip2 lying in H ± . Clearly is a constant times pat_i(0) £jv(x) — £>tv(0) ^at_i(x). Since 
L^ ) (0)/L^1 1 (0) = (N + a)/N we see using ((□]) that we may take 

i/j l = VN£ N - y/N + a£ N -i. 

For -^2; it follows from the discussion near the beginning of the last section that the appro- 
priate linear combination of £jv and £,n-i m ay be found by requiring that it be orthogonal to 
lPq. From the fact that L^\x) x a ~ l e~ x dx = T(a) and from ( |4.1|) we see that the linear 
combination 

does the job. 

Solving for £n and (,n-i in terms of ipi and ip2 and substituting into ( |4~2"|) we obtain for 
the kernel of [D, K], 



^N(N + a) 



ipi(x) if} 2 (x) 
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Therefore 

y/N(N + a) 



X 



Hence for this ensemble we find that S N (x,y) is equal to K^(x 1 y) plus 



2 

( 4 ), 



N{N + a) / ,— l 

- (JnTo^ i N {x) - VN Zn-^x)) (Vn e£ N (y) - VnT^ e^-i(y)) (4.3) 



2 

7(1), 



and that S N (x,y) is equal to K N (x : y) plus 
\JN(N + a 



(VNZ N (x) - VN + a Civ-i(x)) (VN + a e£ N (y) - VN e£ N -i(y)). (4.4) 

These were established for a > 0. For — 1 < a < we must find the analytic con- 
tinuations of the factors in ( f4.3|) and (|4.4|) . The first factors cause no difficulty since they 
are defined and analytic for a > —1. The same is true of the second factor in (fOI) since 



'N pn — \/N + apN-i has zero constant term for a > (since ipi G H) and so for all a. 
The second factor in ( |4.4| ) requires analytic continuation. Assuming at first that a > 
we write it as 

) 1 roc 

(VN + a £ N {z) - VN £at-i(z) dz + - / (VN + a £ N (z) - VN £ N -i(z)) dz. 

2 Jo 

Now the fact eipx G 7i established earlier implies that J °° ipi(z) dz = 0. This is equivalent to 



N i N (z)dz = VN + a / i N ^(z)dz 



so the last integral, with its factor 1/2, is equal to 
Hence the second factor in (|4.4|) is equal to 



/•oo 1 rca 

~ (VN + ^S N (z)-VNZ N -i(z))dz- / z a ' 2 (e- z l 2 VN (z))'dz. (4.5) 

Jy VN + a Jo 

This is analytic for all a > — 1 and so provides the desired analytic continuation. 
When a = 

Ljv(x) - Ljv-i(z) = —L' N (x), 

and we find that 

S%\x,y) = K N (x,y) + V*/ 2 L' N (x) f ' e~ z ' 2 L' N (z) dz, 

I Jo 

S$\x,y) = K N (x,y) + -e^ 2 L' N (x) ( £ e^ 2 L' N (z) dz + l) . 
For both we used the fact J °° ij)\(z) dz = once again and, for the latter, (p~5|). 
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